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摘要 在 不 同 的 轨道 预报 场景 中 , 使 用 的 动力 学 模型 也 不 同 . 例如 , 在 低 轨 空间 碎片 的 预报 中 大 气 阻力 是 十 分 
重要 的 摄 动力 , 而 到 了 中 高 轨 , 大 气 阻 力 就 可 以 忽略 不 计 . 如 何 为 不 同 轨道 类 型 的 空间 碎片 选择 最 优 (满足 精度 
要 求 下 的 最 简 ) 动 力学 模型 还 没有 系统 、 详 尽 的 研究 . 将 对 不 同 精度 需求 、 不 同 轨道 类 型 下 的 大 批量 轨道 进行 预 
报 , 通过 比较 不 同 动力 学 模型 下 的 预报 结果 , 给 出 各 种 预报 场景 的 最 优 动 力学 模型 建议 . 可 以 为 不 同 轨 道 类 型 的 
空间 碎片 在 轨道 预报 时 选择 基准 动力 学 模型 提供 参考 或 标准 . 
关键 词 天 体力 学 , 地 球 , 方法 : 数值 , 行星 和 卫星 : 动力 学 演化 和 稳定 性 
中 图 分 类 号 : P135; ”文献 标识 码 : A 
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1 引言 动力 大 小 , 从 图 中 我 们 可 以 看 到 , 随 着 轨道 高 度 的 
截至 2021 年 ， 太 空中 有 超过 36000 个 尺寸 大 于 变化 , 各 项 摄 动力 的 量 级 有 着 巨大 的 变化 , 起 主导 
10 cm 的 空间 碎片 趾 ， 它 们 包括 废弃 的 卫星 、 火 作用 的 摄 动 力也 在 不 断 变化 . 在 保证 预报 精度 的 前 


日 -二 


箭 残 骸 等 ,而 撞击 事件 、 自 爆 事 件 更 会 使 空间 碎 提 下 , 尽 可 能 地 简化 动力 学 模型 可 以 很 大 地 提高 


) 这 


片 的 数目 爆炸 式 地 增长 . 大 量 聚 集 的 空间 碎片 会 批量 空间 碎片 轨道 预报 的 效率 , 如 何在 保证 精度 和 
威胁 在 轨 航 天 器 的 安全 ， 并 且 严 重 影响 天 文学 观 保证 速度 之 间 取 得 平衡 正 是 这 篇 文章 研究 的 重点 . 
测量 , 太空 中 的 碰撞 事件 已 发 生 多 次 ,比如 2009 年 北美 防空 司令 部 (The North American Aero- 
pT 俄罗斯 的 退役 卫星 宇宙 2231 号 和 美国 的 商业 通讯 space Defense Command, NORAD) 通 过 两 行 根 数 
和 卫星 久 星 33 发 生 了 碰撞 事故 , 二 者 在 碰撞 后 解体 成 (Two-Line Element Set，T 了 TLE) 编目 空间 目标 ， 
为 3000 多 个 空间 碎片 外， 因为 空间 碎片 数量 巨大 ， 并 使 用 SGP4 (Simplified General Perturbations 
分 布 十 分 广泛 , 所 以 在 对 它们 进行 轨道 预报 的 过 4)/SDP4 (Simplified Deep-Space Perturbations 
程 中 , 如 何 选取 适当 的 动力 学 模型 至 关 重 要 , 精细 4) 算 法 对 目标 库 进 行 应 用 . 这 种 算法 是 一 种 基于 正 
的 摄 动 模型 可 以 带 来 精度 更 高 的 结果 , 但 随 之 而 来 则 变换 的 轨道 分 析 解 , 由 于 SGP4 和 SDP4 适 用 的 轨 
的 是 更 大 的 计算 量 需 求 . 对 于 不 同 高 度 、 不 同 偏心 道 不 同 , 它们 考虑 的 动力 学 模型 也 不 同 , 前 者 考虑 
率 的 空间 碎片 , 使 用 的 动力 学 模型 差别 很 大 . 图 1 展 是 地 球 非 球形 引力 位 的 低 阶 近似 以 及 符合 寡 律 
示 了 典型 面 质 比 的 空间 碎片 在 不 同 高 度 受到 的 摄 分 布 的 大 气 阻 力 , 而 后 者 考虑 的 是 地 球 的 非 球形 引 
2022-04-18 收 到 原稿 , 2022-06-06 收 到 修改 稿 
* 国 家 自然 科学 基金 项 目 (11873031) 资 助 


ijstang@nju.edu.cn 


| 


吧 


lml 


到 辣 


64 卷 天 


力 位 的 低 阶 近 似 和 日 月 引力 回 . 然而 这 种 算法 使 用 
的 动力 学 模型 仅 能 保证 特定 的 精度 (在 一 周 内 的 轨 
道 预报 可 以 保证 位 置 误差 在 10 km 量 级 引 5)， 并 且 
依赖 于 其 发 布 的 TLEI. 针对 空间 碎片 , 需要 提出 
一 种 针对 不 同 预报 场景 ( 即 不 同 轨道 类 型 ,不 同 精 
度 需 求 ) 的 动力 学 模型 . 


Acceleration/(m-.s’?) 


Solar Radiation Pressure 


\ 
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图 1 不 同 高 度 的 空间 碎片 受到 的 摄 动力 . 其 中 GM 表示 地 球 的 二 体 引 
力 , J 或 J 表示 地 球 的 非 球形 引力 摄 动 的 主要 带 谐 项 和 主要 田 谐 项 ， 
计算 大 气 阻力 和 光 压 摄 动 的 相关 参数 见 2.1 节 . 


Fig.1 The perturbation force of space debris at different 
heights. Among them, GM represents the two-body 
gravitation of the earth, J, or J,: represents the main band 
and tesseral harmonic terms of the aspheric gravitational 
perturbation of the earth. The relevant parameters for 
calculating the atmospheric drag and solar radiation 


pressure perturbation are listed in Section 2.1. 


我 们 使 用 不 同 的 动力 学 模型 对 大 量 的 空间 碎 
片 轨道 进行 轨道 预报 , 通过 比较 不 同 模型 下 轨道 预 
报 的 结果 ， 以 精细 动力 学 模型 下 的 结果 为 参考 ,给 
出 不 同 的 动力 学 模型 下 轨道 预报 的 精度 . 参照 预报 
结果 的 精度 , 我 们 可 以 针对 不 同 的 预报 场景 , 给 出 
最 优化 的 动力 学 模型 . 

本 文 第 2 节 将 会 描述 轨道 预报 的 过 程 , 第 3 节 会 
展示 不 同 预 报 场景 和 不 同 动力 学 模型 下 的 预报 精 
度 , 第 4 节 通 过 分 析 第 3 节 的 结果 给 出 动力 学 模型 建 
议 , 第 5 节 是 总 结 . 
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2 ”轨道 预报 方法 
2.1 ”空间 碎片 的 动力 学 

空间 碎片 在 地 球 的 引力 下 运动 , 在 理想 情况 下 
为 一 个 二 体 运 动 , 它 是 一 个 可 积 系统 , 存在 解析 解 . 
而 在 现实 中 , 地 球 质量 和 形状 的 不 均匀 、 太 空中 稀 
薄 大 气 的 存在 等 因素 使 得 空间 碎片 绕 地 球 运动 成 
为 一 个 受 摄 动 的 二 体 问 题 . 我 们 研究 的 摄 动力 主要 
包括 地 球 的 非 球形 引力 、 大 气 阻 力 、 日 月 的 第 三 
体 引力 、 太 阳 辐 射 压 力 (Solar Radiation Pressure， 
SRP)、 地 球 固体 潮汐 力 、 后 牛顿 效应 和 太阳 系 其 
他 大 天 体 的 第 三 体 引力 . 

本 文中 使 用 的 地 球 引 力 场 模型 是 GGM03C 
(Global Gravity Model 03C)W， 它 是 由 GRACE 


(Gravity Recovery and Climate Experiment) 卫 星 


数据 反 演 的 引力 场 之 一 ， 最 高 支持 360 x 360 阶 次 . 
在 计算 第 三 体 引力 摄 动 时 , 计算 日 月 位 置 使 用 的 是 
JPL (Jet Propulsion Laboratory) 数 值 历 表 和 一 种 
平均 化 的 分 析 历 表 , 我 们 还 会 在 下 文中 对 二 者 做 一 
个 简单 的 比较 . 在 大 气 阻 力 摄 动 的 计算 中 , 相对 麻 
烦 的 是 大 气 密度 的 计算 , 因为 大 气 受 太阳 辐射 强度 
和 地 磁 活 动 等 影响 较 大 ， 大气 密 度 模 型 的 精度 有 
限 由 ,本 文中 使 用 的 是 一 种 参考 了 NRLMSISE-00 
(Naval Research Laboratory Mass Spectrometer 
and Incoherent Scatter Radar Exosphere 2000)I9 
大 气 模型 的 扁 球 指数 平均 大 气 模型 央 . 大 气 阻力 和 
光 压 摄 动 的 强度 受 面 质 比 等 因素 的 影响 , 在 没有 特 
别 说 明 的 情况 下 , 本 文采 用 表 1 所 示 参 数 . 表 1 中 列 
举 的 4 项 参数 分 别 为 大 气 阻力 系数 、 迎 风向 等 效 面 
质 比 、 太 阳 方 向 等 效 面 质 比 和 反射 系数 , 它们 将 被 
应 用 到 大 气 阻 力 摄 动 和 太阳 辐射 压 摄 动 的 计算 中 . 

轨道 预报 采用 数值 方法 ， 即 使 用 变 步 长 的 
Runge-Kutta-Felhberg7 (8) 积 分 器 积分 高 斯 型 摄 动 
运动 方程 , 截断 误差 小 于 10-93 (长 度 单位 为 地 球 半 
径 , 角度 单位 为 弧度 ), 轨道 根 数 采用 第 1 类 无 奇 点 
根 数 , 妈 a,i, ,EE = ecosw,n = 二 esinw, 和 =w+M, 
其 中 o 为 轨道 半 长 径 、? 为 轨道 倾角 、9 为 升 交 点 经 
度 、e 为 轨道 偏心 率 、w 为 近地点 辐 角 、M 为 平 近 
点 角 . 
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表 1 大 气 阻力 和 太阳 辐射 压 相 关 人 参数 取 值 表 
Table 1 Table of parameters related to 


atmospheric resistance and solar radiation 


pressure 


Parameters Values 


Drag coefficient 2.2 


Windward effective et 
7.08x10 ~ m :kg 

surface-to-mass ratio 
Effective surface-to-mass 2 _1 
7.08x10 ~ m’ .kg 
ratio towards the Sun 


Reflection coefficient 0.5 


2.2 ”试验 轨道 的 初 值 选择 

为 了 模拟 真实 的 空间 碎片 分 布 , 轨道 初 值 将 会 
根据 TLE 数 据 进 行 选 择 , 一 方面 使 其 符合 真实 的 空 
间 碎 片 分 布 (主要 是 半 长 径 和 偏心 率 的 分 布 ), 另 一 
方面 需 保证 均匀 且 全 面 地 取 值 如 图 2 所 示 , 我 们 
画 出 了 目前 TLB 所 有 在 编 空间 目标 的 轨道 高 度 与 
轨道 偏心 率 分 布 , 并 且 使 用 4 种 不 同 的 颜色 代表 上 
述 4 种 轨道 类 型 . 我 们 对 空间 目标 的 轨道 进行 划分 ， 
分 为 以 下 4 类 ,这 里 轨道 高 度 的 定义 为 轨道 的 半 


oo Typel 
人 Type2 
D Type3 
9 Type4 


Eccentricity 


0 
102 103 


长 径 与 地 球 赤 道 半径 之 差 , 即 轨道 近地点 与 远地点 
的 平均 高 度 . 

类 型 一 : 低 轨 (Low Earth Orbit, LEO), 轨道 
高 度 在 2000 km 以 下 , 其 中 比较 具有 代表 性 的 是 太 
阳 同 步 轨道 和 星 链 卫星 所 在 轨道 . 我 们 研究 的 低 轨 
目标 总 数 为 15921 个 . 

类 型 二 : 中 轨 (Medium Earth Orbit, MEO) 小 
偏心 率 , 轨道 高 度 高 于 2000 km, 低 于 36000 km，e 
< 0.1. 我 们 研究 的 中 轨 小 偏心 率 目标 总 数 为 584 个 . 

类 型 三 : 中 轨 大 偏心 率 ， 轨 道 高 度 高 于 
2000 km， 低 于 36000 km，e> 0.1,， 包括 闪电 轨道 
(Molniya). 我 们 研究 的 中 轨 大 偏心 率 目 标 总 数 为 
2264 个 . 

类 型 四 : 高 轨 (High Earth Orbit， HEO)， 
一 部 分 为 地 球 同 步 轨道 (Geosynchronous Or- 
bit，GSO), 其 中 包括 倾斜 地 球 同步 轨道 (Inclined 
Geosynchronous Orbit， IGSO)、 地 球 静 止 轨 
道 (Geostationary Orbit，GEO)， 另 一 部 分 为 高 度 
高 于 GSO 的 轨道 (包含 大 偏心 率 轨道 和 小 偏心 率 轨 
道 ) 和 高 度 等 于 GSO 的 大 偏心 率 轨道 . 我 们 研究 的 
高 轨 目 标 总 数 为 1560 个 . 


Height/km 


图 2 ” TLE 在 编目 标 随 轨道 高 度 与 轨道 偏心 率 的 分 布 


Fig.2 The distribution of TLE objects with orbit height and orbit eccentricity 
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3 ”轨道 预报 的 误差 分 析 
在 以 下 小 节 中 , 我 们 比较 的 都 是 一 天 内 的 位 置 


偏差 . 假设 待 研究 的 动力 学 模型 为 A, 基准 动力 学 模 
型 为 B. 基准 动力 学 模型 中 的 摄 动力 计算 是 足够 精 
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一 个 处 于 中 轨 的 小 偏心 率 轨道 的 空间 碎片 , 那么 引 
力 场 模型 就 可 以 简化 到 10x10 阶 , 对 于 地 球 同 步 轨 
道 的 空间 碎片 , 更 是 可 以 简化 到 5x5 阶 . 

对 于 低 轨 空间 碎片 的 轨道 预报 ， 引 力 场 的 准 
确 性 要 求 较 高 , 随 着 引力 场 阶 次 的 提高 , 轨道 预报 
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的, 壁 如 在 研究 10 x 10 非 球形 引力 位 对 预报 精度 
的 影响 时 , 动力 学 模型 A 中 的 引力 场 阶 次 为 10 x 10， 
而 基准 动力 学 模型 B 中 的 引力 场 阶 次 为 80 x 80. 分 
别 使 用 这 两 种 动力 学 模型 预报 同一 条 轨道 ， 每 隔 

一 分 钟 记录 一 次 , 得 到 两 条 星 历 . 用 位 置 表示 两 条 
星 历 分 别 为 7 ， TA2, 人 
TB1440: 定义 位 置 偏差 4 为 


/> |r — Fil? 
1440 

这 个 偏差 就 表征 了 使 用 动力 学 模型 A 在 一 
的 位 置 误差 . 
引力 场 阶 次 对 轨道 预报 精度 的 影响 
作为 地 球形 状 和 质量 分 布 不 均匀 导致 的 摄 动 
力 , 非 球形 引力 摄 动 的 效应 会 随 轨道 高 度 的 增加 而 
减弱 . 另外 , 考虑 的 引力 场 阶 次 越 高 , 摄 动力 的 计算 
越 精 确 , 但 是 轨道 预报 的 速度 也 越 慢 . 不 同 轨道 类 
型 对 引力 场 的 阶 次 的 敏感 程度 不 同 , 为 了 比较 不 同 
引力 场 阶 次 下 的 位 置 偏差 , 我 们 使 用 不 同 精 细 程 度 
的 引力 场 对 空间 碎片 进行 轨道 预报 , 除了 地 球 的 非 
球形 引力 摄 动 外 , 这 里 没有 考虑 别 的 摄 动 力 . 用 以 
计算 的 初始 轨道 是 2.2 节 中 提 到 的 全 部 轨道 . 

不 同 引 力 场 阶 次 下 的 轨道 预报 误差 如 图 3-4 所 
示 ， 其 中 图 3 对 应 低 轨 目标 , 图 中 每 一 个 点 表示 相 
应 轨道 高 度 和 引力 场 阶 次 下 的 位 置 误差 位 置 误 
差 的 大 小 用 颜色 表示 . 线条 为 特定 精度 的 分 界线 ， 
Fe os li te tld, 分 界 
线 上 方 区 域 的 精度 高 于 10 m， 分 界线 下 方 区 域 的 
精度 低 于 10 m. 图 4 对 应 中 高 轨 目 标 , 图 中 的 表示 
方法 与 图 3 类 似 . 从 图 3-4 中 可 以 看 出 , 地 球 非 球形 
引力 摄 动 的 影响 随 着 轨道 高 度 的 增加 而 降低 ， 这 
就 意味 着 预报 的 轨道 越 高 ,采用 引力 场 模 型 就 可 
以 越 简 单 . 例如 , 我 们 希望 轨道 预报 的 结果 在 一 天 
内 的 位 置 误差 不 超过 100 m， 如 果 预 报 的 对 象 是 一 
个 高 度 为 300 km 的 低 轨 空 
330x30B 
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) TA1440 和 和 7B1， TB2, "**， 


(1) 
天 内 积累 


3.1 
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介 的 引力 场 . 相同 精度 要 求 下 ， 如 果 这 


42-4 


位 置 偏差 的 减 小 十 分 缓慢 ， 当 引力 场 采用 80x80 阶 
时 , 绝 大 部 分 空间 碎片 的 轨道 预报 位 置 偏差 减 小 到 
了 1m 以 下 , 但 这 样 高 阶 次 引力 场 的 计算 量 已 经 十 
分 大 了 , 再 提高 轨道 预报 的 精度 十 分 困难 . 当 轨 道 
较 低 (例如 近地点 低 于 400 km 的 轨道 ) 时 , 随 着 引力 
场 阶 次 的 提高 , 精度 提升 十 分 缓慢 . 如 果 我 们 希望 
这 些 空间 碎片 的 轨道 预报 达到 厘米 级 精度 , 这 将 会 
是 十 分 困难 的 . 另外 , 这 样 的 精度 要 求 也 并 不 现实 ， 
因为 在 低 轨 空间 碎片 的 轨道 预报 中 , 大 气 阻 力 摄 动 
的 影响 显著 ， 且 精确 计算 十 分 困难 , 所 以 由 大 气 模 
型 造成 的 误差 一 般 远 大 于 1 cm. 因此 , 对 于 轨道 高 
度 低 于 400 km 的 轨道 , 将 引力 场 的 精度 定 在 厘米 或 
以 下 意义 不 大 . 事实 上 , 大 多 数 低 轨 空间 碎片 不 会 
处 在 极 低 的 轨道 上 , 随 着 轨道 高 度 的 增加 , 对 引力 
场 精确 度 的 依赖 逐渐 在 减 小 . 同样 很 难 使 精度 提高 
到 厘米 级 的 是 中 轨 大 偏心 率 轨道 , 其 中 一 些 轨 道 的 
近地点 高 度 甚至 低 于 200 km, 所 以 对 他 们 的 处 理 与 
低 轨 类 似 . 
使 用 不 同 的 引力 场 阶 次 需要 的 算 力 有 很 大 区 
别 , 因此 选取 合适 的 引力 场 是 节省 轨道 预报 时 间 的 
关键 . 如 图 5 所 示 我 们 展示 了 使 用 不 同 精细 程度 的 
引力 场 预报 一 天 轨道 所 需 的 时 间 . 使 用 不 同 尼 
或 软件 环境 时 ,所 需 时 间 会 有 所 差异 , 我 们 这 里 
用 的 CPU 是 intel core i5-9500, 使 用 的 编译 器 是 in- 
tel visual fortran, 使 用 release 方 式 编 译 . 图 中 展示 
的 4 条 轨道 的 初始 开 普 勒 根 数 如 表 2 所 示 . 图 5 中 纵 
坐标 表示 轨道 预报 一 天 使 用 的 CPU 时 间 , 取 10 次 测 
试 的 平均 值 . 一 般 来 说 轨道 越 高 , 轨道 预报 的 速度 
越 快 , 这 是 因为 轨道 高 度 越 高 , 变 步 长 积分 器 可 以 
取 到 的 步 长 越 长 , 同样 的 预报 时 间 内 步 数 可 以 越 少 . 
另 一 方面 , 对 于 同一 高 度 的 轨道 来 说 , 动力 学 模型 
越 简单 , 轨道 预报 速度 越 快 . 从 图 5 中 的 LEO 轨 道 来 
看 , 如 果 将 引力 场 阶 次 从 10x10 提 升 到 50x50, 轨道 
预报 的 用 时 增加 到 了 原来 的 5 倍 左右 . 因此 , 如 果 结 
合 图 3 与 图 4 给 出 某 种 预报 场景 下 的 最 佳 引力 场 阶 
次 , 就 可 以 大 大 地 节省 算 力 . 
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图 3 低 轨 空间 碎片 在 不 同 阶 次 的 引力 场 下 预报 结果 的 位 置 误 差 . 横 坐 标 为 轨道 的 平均 高 度 ， 纵 坐标 为 引力 场 的 阶 次 , 其 中 2X0 表 示 仅 考虑 .7 项 , 线 
条 为 某 一 个 位 置 误差 的 分 界线 . 


Fig.3 The position error of the LEO space debris projected in different orders of gravitational field. The abscissa is the average 
height of the orbit, and the ordinate is the order of the gravitational field, where 2x0 means that only .12 is considered, and the 


line is the dividing line of a certain position error. 
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图 4 中 高 轨 空 间 碎 片 在 不 同 阶 次 的 引力 场 下 预报 结果 的 位 置 误差 . 横 坐 标 为 轨道 的 平均 高 度 , 纵 坐 标 为 引力 场 的 阶 次 , 其 中 2x0 表 示 仅 考虑 .7 项 ， 
线条 为 某 一 个 位 置 误差 的 分 界线 . 


Fig.4 The position error of the projection results of MEO/HEO space debris under different orders of gravitational field. The 


abscissa is the average height of the orbit, and the ordinate is the order of the gravitational field, where 2x0 means that only .72 


is considered, and the line is the dividing line of a certain position error. 
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到 5 轨道 预报 用 时 随 引力 场 阶 次 的 变化 图 . 当 
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当 横 坐标 为 0 时 ,表示 二 体 引力 ， 
闪电 轨道 是 一 种 特殊 的 大 偏心 率 轨道 , 轨道 周期 约 为 12 h. 


当 横 坐标 为 n 时 , 表示 n xn 阶 次 引力 场 . 图 中 Molniya 表 示 闪 电 轨 道 ， 


Fig.5 The orbit prediction time as a function of the order and the degree of the gravitational field. When degree or order is 0， 


it represents two body gravity. When degree or order is n, it represents the gravitational field of order nxn. Molniya in the 


figure represents Molniya orbit. Molniya orbit is a special orbit with large eccentricity, and the orbital period is about 12 h. 


表 2 4 条 轨道 在 初始 时 刻 的 开 普 勒 根 数 


Table 2 Keplerian elements of the four orbits at 


the initial moment 


Orbit a/km e 1 Qf w/ MA 
LEO 7178 0.01 45 30 40 50 
MEO 10378 0.01 45 30 40 50 
Molniya 26600 0.74 63 30 270 50 
GEO 42378 0.01 0 30 40 50 
3.2 ”大 气 阻力 对 轨道 预报 精度 的 影响 
本 节 对 低 轨 试验 轨道 进行 一 天 的 轨道 预报 , 对 
是 否 考虑 大 气 阻力 的 预报 结果 , 研究 动力 学 模型 


A ee lo 
. 现 有 研究 表明 , 目 


中 , 大气 密 度 的 计算 尤其 重要 


前 的 大 气 密 度 模 型 在 比较 算 例 中 各 有 从 


来 计算 的 轨道 是 


劣 ， 
之 间 内 符合 误差 最 大 可 达 20%-30%00-12]. ] 
类 型 一 的 轨道 , 即 低地 球 轨道 , 总 


模型 


明 的 情况 下 , 我 们 取 太 阳 活 动 指数 F107 = 130, 地 
球 磁场 指数 Ap==10. 在 这 种 环境 下 , 高 度 440 km 的 
大 气 密度 约 为 1.627x10-12 kg.m-3， 在 本 小 节 我 
们 也 会 适当 研究 不 同 F107 和 Ap 下 的 位 置 偏差 的 区 


别 . 


轨道 预报 的 结果 如 图 6 所 示 , 这 里 横 坐 标 取 近 
地 点 高 度 , 这 是 因为 即使 是 半 长 径 很 大 的 轨道 , 如 
人 
气 , 使 得 轨道 受到 大 气 的 影响 依旧 十 分 明显 , 所 以 
使 用 近地点 高 度 衡量 大 气 阻力 的 地 位 更 为 合理 . 黑 
色散 点 表示 在 一 般 情况 下 (F107 = 130, Ap = 10) 
试验 轨道 的 误差 情况 ， 实 线 是 对 黑色 散 点 的 拟 合 
这 里 使 用 的 拟 合 公式 是 y = e- 订 六 ,使 用 这 种 公式 
构 型 的 原因 是 这 里 的 位 置 偏差 与 大 气 密度 存在 密 
切 的 关系 ， 而 大 气 密 度 随 高 度 的 变化 正 是 这 种 形 


这 里 


数 1600 条 . 随 着 太阳 活动 强度 、 地 球 磁 场 等 环境 的 
变化 , 大 气 密度 的 改变 也 相当 显著 , 在 没有 特别 说 


式 , 即 一 方面 大 气 密度 随 高 度 指数 式 下 降 , 另 一 方 
面 下 降 的 速度 会 越 来 越 慢 (标高 越 来 越 高 ). 在 一 天 
的 轨道 预报 中 ，400 km 左右 高 度 的 空间 碎片 可 以 
职 累 几 公里 到 几 十 公里 的 误差 ， 而 当 轨 道 高 度 低 
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于 400 km 时 ， 大 和气 阻力 在 一 天 内 对 预报 结果 的 影 
响 会 呈现 指数 式 的 上 升 ， 甚 至 可 以 达到 百 公里 和 


干 公 里 的 量 级 . 同时 , 这 部 分 轨道 的 高 度 下 降 也 是 
十 分 快 的 , 因此 空间 碎片 或 航天 器 所 处 的 轨道 高 度 


一 般 不 会 低 于 300 km. 而 随 着 轨道 稍稍 抬 高 , 大 气 
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阻力 的 影响 开始 又 减 . 几 百 公里 高 的 轨道 积累 的 
误差 只 有 几米 到 几 十 米 之 间 ， 对 于 轨道 高 度 大 于 
1200 km 的 空间 碎片 , 位 置 误差 均 小 于 1 m 且 迅速 误 
减 ， 到 了 1600 km， 则 一 般 场 景 下 几乎 可 以 不 用 考 


一 = Ap = 5, F107 = 70 


uaan Ap= 15,F107= 190 
© SSO or StarLink 


cs 六 入 [moo 
| 中 
二 
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图 6 大 气 阻 力 对 轨道 预报 精度 的 影响 随 轨道 高 度 的 变化 图 . 方形 散 点 为 相 


应 情况 下 不 同 轨道 的 位 置 偏差 的 分 布 情况 , 实 线 是 对 方形 散 点 的 拟 合 . 虚 


线 和 点 划 线 是 其 他 两 种 情况 下 的 结果 的 拟 合 . 圆 形 散 点 表示 的 是 低 轨 


标 


bh 的 一 些 具 有 代表 性 的 轨道 , 太阳 同步 轨道 (Sun-Synchronous Orbit， 


SSO) 和 星 链 (StarLink) 卫 星 轨道 . 


Fig.6 The effect of atmospheric resistence on the accuracy of orbit prediction as a function of orbital altitude. The square 


scatters are the distribution of the positional deviations of different orbits in the corresponding cases, and the solid line is the 


fitting of the square scatter points. The dashed and dot-dashed lines are the fittings of the results in the other two cases. The 


circle scatters represent SSO and starlink satellite orbits. 


SSO 卫 星 和 星 链 (StarLink) 卫 星 等 是 在 低 轨 目 
标 中 的 一 些 典 型 目标 , 我 们 在 图 6 中 标记 了 这 些 目 
标 . 从 图 中 我 们 可 以 看 出 , 这 些 目 标 相 比 其 他 目标 
并 没有 特殊 之 处 , 即 决 定 大 气 阻 力 造 成 的 误差 情况 
主要 在 于 轨道 的 平均 高 度 或 近地点 高 度 . 

我 们 还 在 图 6 中 展示 了 其 他 两 种 太阳 活动 情况 
和 地 磁 活 动情 况 下 大 气 阻力 对 轨道 预报 结果 的 影 
响 情 况 , 这 两 种 情况 分 别 是 F107 = 70，Ap = 5 和 
下 107 = 190，Ap = 15. 因为 过 多 的 散 点 会 重 对 在 
一 起 使 图 像 过 于 混乱 , 这 里 没有 画 出 这 两 种 F107 和 
Ap 下 相应 的 散 点 图 而 仅 展 示 了 拟 合 结果 , 拟 合 结 
果 分 别 用 点 划 线 与 虚线 表示 . 拟 合 公式 同样 是 取 
y 二 e- 专 专 的 形式 , 拟 合 参数 如 表 3 所 示 . 太阳 活动 


和 地 磁 活 动 会 影响 大 气 密度 , 从 而 间接 地 影响 大 气 
阻力 的 大 小 , 并 且 影 响 程度 比较 显著 . 就 轨道 高 度 
为 800 km 的 空间 碎片 而 言 , 虚线 与 点 划 线 的 结果 有 
一 个 量 级 左右 的 差别 . 因此 , 在 考虑 大 气 阻 力 对 轨 
道 预 报 结果 的 影响 时 , 有 必要 考虑 当日 的 太阳 活动 
情况 和 地 磁 活 动情 况 . 

3.3 ”日 月 第 三 体 引 力 对 轨道 预报 精度 的 影响 

第 三 体 引 力 摄 动 ( 这 里 指 太阳 和 月 球 ) 是 空间 碎 
片 动力 学 中 的 主要 摄 动力 来 源 之 一 ， 在 绝 大 多 数 
的 预报 场景 中 都 不 可 忽略 , 对 于 任意 高 度 的 空间 
碎片 , 太阳 和 月 球 的 引力 不 仅 在 同一 个 量 级 , 而 且 
大 小 接近 , 那么 我 们 就 不 需要 对 他 们 分 别 研 究 , 而 


| 
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是 考虑 它们 的 共同 作用 . 第 三 体 引力 摄 动 的 影响 会 预报 场景 的 需求 . 

随 着 空间 碎片 轨道 高 度 的 增加 而 增加 , 对 于 低 轨 空 

间 碎 片 而 言 ， 合 不 考虑 第 三 体 引 力 摄 动 的 动力 学 Table 3 dt epi the three 
模型 预报 一 天 的 弧 段 会 产生 百 米 级 的 位 置 误差, 而 ee 

对 于 地 球 同步 轨道 空间 碎片 而 言 , 位 置 误差 会 达到 nn. 

几 十 公里 . 第 三 体 引 力 摄 动 的 计算 并 不 复杂 , 关键 


130 10 7404 一 4.323 201.4 


在 于 获得 准确 的 太阳 、 月 球 位 置 . 我 们 将 会 使 用 
两 种 方法 来 获得 太阳 和 月 球 位 置 , 分 别 是 分 析 历 


70 5 4678 ”一 3.050 48.12 


a 表 鸣 ). 在 这 一 小 节 我 们 190 15 10930 —5.817 423.8 
要 研究 两 个 问题 , 其 一 是 分 析 历 表 和 数值 历 表 的 
精度 对 轨道 预报 的 时 光 啊 有 多 大 ? 其 二 是 如 果 使 用 图 7 展示 了 使 用 日 月 分 析 历 表 带 来 的 位 置 误差 ， 


不 同 版 本 的 JPL 历 表 , 它们 的 差别 有 多 大 ? 分析 历 即 相 对 于 使 用 精度 十 分 高 的 数值 历 表 结 果 的 位 置 
表 中 的 太阳 位 置 可 以 使 用 太阳 轨道 的 平 根 数 获得 ， 偏差 . 这 里 计算 的 轨道 分 为 4 部 分 , 分 别 是 低地 球 
而 月 球 位 置 的 周期 项 较 大 , 不 可 以 简单 地 使 用 平 根 轨道 、 中 轨 小 偏心 率 轨 道 、 闪 电 轨 道 和 地 球 同步 
数 获 得 , 在 本 文中 使 用 的 是 带 有 部 分 周期 项 的 解析 轨道 . 显然 , 轨道 高 度 越 高 的 轨道 受到 第 三 体 引 力 
轨道 1 站, 相对 精度 (即位 置 误差 与 位 置 大 小 之 比 ) 约 摄 动 的 影响 越 大 , 分 析 历 表 带 来 的 误差 也 越 大 . 对 
为 |Ar|/r' = O(10-3). 分 析 历 表 不 仅 计算 速度 快 ， 于 处 于 低 轨 的 空间 碎片 , 一 天 内 积累 的 位 置 误差 
而 且 不 依赖 于 大 体 量 的 拟 合 参数 数据 文件 , 可 以 减 在 1 m 左 右 , 对 于 处 于 中 轨 的 空间 碎片 , 一 天 内 积累 
小 程序 的 内 存 占 用 , 提高 运行 速度 . 相对 于 JPL 历 的 位 置 误差 在 10 m 量 级 , 对 于 闪电 轨道 和 地 球 同步 
表 , 分 析 历 表 精 度 较 低 , 但 依然 能 满足 大 多 数 轨道 轨道 的 空间 碎片 , 位 置 误差 则 在 百 米 量 级 . 
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妈 7 ”使 用 日 月 分 析 历 表 产 生 的 位 置 误差 随 轨道 类 型 的 变化 图 . 从 左 到 右 按 轨道 高 度 排序 , 依次 是 低 轨 、 中 轨 近 圆 、 闪 电 轨 道 和 地 球 同步 轨道 . 


Fig.7 The position error caused by the analytical ephemeris of the sun and the moon is plotted against the type of orbit. From 
left to right, they are sorted by orbital height, followed by low orbit, medium orbit near circular, Molniya orbit, and 


geosynchronous orbit. 
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JPL 历 表 是 由 美国 喷气 推进 实验 室 研 究 并 发 布 置 精 度 ， 以 最 新 的 DE440 为 基准 . 图 8 ( 左 ) 为 月 球 
的 基于 多 项 式 拟 合 的 行星 历 表 . 最 早 获 得 广泛 应 用 的 位 置 精度 ， 图 8 ( 右 ) 表 示 太 阳 的 位 置 精度 ， 纵 
的 是 DE96H51， 后 续 JPL 历 表 发 展 了 许多 版 本 ， 比 坐标 为 日 月 位 置 的 相对 误差 ， 即 e = |7bExw 一 
如 1985 年 的 DE125 根 据 旅行 者 二 号 飞 掠 天 王 星 的 7bp440|/|7bpaso|， 其 中 7zDpywx 是 相应 历 表 下 的 日 月 
数据 对 历 表 精度 进行 了 修正 19. 本 文中 轨道 预报 默 位 置 ， 例 如 Fpaso，7Fibr200， 月 球 的 相对 位 置 误差 
认 使 用 的 是 DE430H4, 这 是 目前 比较 稳定 且 精 度 比 从 DE200 的 10-6 降 低 到 了 到 DE430 的 10-3. 太阳 的 
较 高 的 版 本 , 而 最 新 的 JPL 历 表 是 DE440. 相对 位 置 误差 从 DE200 的 10-" 降 低 到 了 DE430 的 

如 图 8 展示 的 几 个 典型 版 本 JPL 历 表 的 日 月 位 10-9, 这 也 得 益 于 月 球 激光 测 距 的 发 展 071. 
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图 8 主要 版 本 的 JPL 历 表 相 对 DE440 (Development Ephemeris 440) 的 位 置 偏差 示意 图 


Fig.8 Schematic diagram of the relative position deviation of the major version of the JPL planetary ephemerides compared 
with DE440 


日 月 的 位 置 误差 间接 导致 的 轨道 预报 位 置 误 这 些 轨道 从 TLE 编 目 中 随机 选取 , 包含 了 2.2 节 所 述 
全 差 如 图 9 所 示 . 空间 碎片 的 初始 轨道 是 一 条 高 度 为 的 各 个 类 型 的 轨道 . 图 10 画 出 了 4 条 有 具有 代表 性 轨 
-= 36000 km 的 近 圆 轨道 , 在 半 个 月 的 轨道 预报 中 , 位 道 的 预报 误差 情况 , 它们 分 别 是 低 轨 空间 碎片 、 中 
bd 置 误差 与 预报 时 间 有 接近 线性 的 关系 . 如 果 采 用 轨 空 间 碎 片 、 闪 电 轨 道 空 间 碎 片 和 地 球 同步 轨道 


Wk 


历 表 DE200， 每 天 积累 的 误差 为 厘米 级 ， 而 使 用 空间 碎片 , 初始 轨道 开 普 勒 根 数 如 表 2 所 示 . 


DE430, 每 天 积累 的 位 置 误差 不 到 1 mm. 结果 表明 ， SRP 导 致 的 位 置 误差 有 明显 的 周期 特性 , 周期 

即使 不 使 用 比较 新 版 本 的 JPL 历 表 , 轨道 预报 的 精 。” 正 十 轨道 周期 . 根据 预报 的 结果 , 我 们 发 现 虽然 太 

度 也 不 会 受到 明显 影响 . 昌 辐 射 压力 的 大 小 随 轨 道 高 度 的 变化 并 不 显著 , 但 

是 SRP 对 高 轨 空 间 碎片 的 摄 动 效果 依然 远大 于 低 

3.4 ”太阳 辐射 压 (SRP) 对 轨道 预报 精度 的 影响 轨 空 间 碎片 . 对 于 中 低 轨 空间 碎片 , 不 考虑 SRP 的 

本 文采 用 的 典型 面 质 比如 表 1 所 示 . 这 里 使 用 。 ”预报 结果 在 一 天 内 的 位 置 误差 约 为 10 m, 对 于 地 球 
不 同 的 动力 学 模型 预报 了 800 条 轨道 来 进行 比较 ， ”同步 轨道 和 闪电 轨道 , 误差 约 为 100 m. 


= 
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妈 9 ”不 同 版 本 的 JPL 历 表 对 轨道 预报 精度 的 影响 


Fig.9 The influence of different versions of JPL ephemeris on the accuracy of orbit prediction 
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图 10 ”太阳 辐射 压 对 4 条 轨道 预报 精度 的 影响 


Fig.10 The influence of solar radiation pressure on the accuracy of the four orbits 


在 计算 SRP 摄 动力 时 ， 地 影 的 处 理 十 分 关键 ， 影 中 , SRP 摄 动 不 起 作用 , 在 半 影 中 , 太阳 被 地 球 部 
当空 间 碎片 处 于 地 影 之 内 时 , 不 会 受到 SRP 的 影响 . 分 遮挡 , 由 未 被 遮挡 的 部 分 提供 光 压 摄 动力 , 即 光 
一 般 情况 下 , 为 了 简化 计算 过 程 , 我 们 将 太阳 考虑 压 摄 动力 的 大 小 由 未 被 遮挡 部 分 的 太阳 视 圆 面 所 
为 一 个 足够 远 且 足够 小 的 点 ,那么 地 影 呈 圆柱 形 . 占 整体 视 圆 面 的 比例 决定 Da. 
在 以 上 的 计算 中 , 地 影 模型 均 是 如 此 处 理 . 而 现实 这 里 将 会 研究 不 同 地 影 模型 的 使 用 会 对 轨道 
中 , 因为 太阳 与 地 球 的 距离 为 有 限 值 , 地 影 的 形状 预报 结果 产生 多 大 的 影响 . 因为 地 影 主要 影响 低 
应 为 圆锥 形 . 另外 , 考虑 到 太阳 并 非 一 个 点 而 是 一 轨 空 间 碎片 的 轨道 而 对 中 高 轨 的 空间 碎片 的 轨道 
个 圆 球 , 地 影 会 被 分 成 本 影 和 半 影 两 个 部 分 . 在 本 ”影响 很 小 ,我们 这 里 以 低 轨 空 间 碎 片 作 为 研究 目 
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标 . 目标 轨道 的 轨道 高 度 分 别 为 400 km (LEO1)、 

700 km (LEO2)、 1000 km (LEO3)、 1300 km 
(LEO4), 偏心 率 取 0.01, 轨道 倾角 均 取 45°, 升 交点 
经 度 、 近 地 点 角 距 和 平 近 点 角 均 分 别 取 30*、40°? 和 
50°. 另外 , 由 于 空间 碎片 进出 地 影 的 时 间 很 短 , 通 
常 只 有 几 秒 钟 , 所 以 为 了 体现 出 进出 地 影 边 缘 时 不 
同 地 影 模 型 的 差别 , 在 数值 方法 轨道 预报 中 需要 将 
步 长 调整 到 很 小 的 水 平 , 我 们 这 里 将 步 长 调整 为 固 
定 的 1 s. 将 步 长 固定 为 1 s 后 会 显著 地 影响 轨道 预 
报 的 速度 , 举例 来 说 , 一 条 700 km 高 度 的 LEO 轨 道 ， 
取 动 力学 模型 为 30x30 阶 次 引力 场 、 考 虑 大 气 阻 力 
和 日 月 第 三 体 引力 , 在 变 步 长 情况 下 将 其 轨道 向 前 
预报 一 天 需要 1.8 s 左 右 (不 同 的 硬件 或 软件 条 件 下 
会 有 一 些 差 异 )， 而 相同 条 件 下 的 固定 步 长 1 s 的 预 
报 需 要 68 s 左 右 , 时 间 消 耗 大 大 增加 . 幸运 的 是 , 从 


Deviation/m 
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我 们 的 结果 来 看 , 绝 大 多 数 预报 场景 不 需要 考虑 如 
此 精细 的 地 影 模型 差别 , 也 就 不 必 如 此 大 程度 地 牺 
牲 预 报 速度 . 

图 11 展 示 了 使 用 柱 形 地 影 与 锥 形 地 影 两 种 情 
况 下 的 轨道 预报 结果 的 差别 ， 纵 坐标 是 使 用 两 种 
地 影 模 型 的 预报 结果 的 位 置 差异 , 或 者 说 以 锥 形 
地 影 为 基准 模型 , 圆柱 形 地 影 产 生 的 预报 误差 . 我 
门 可 以 看 到 , 随 着 时 间 的 推移 , 使 用 不 同 的 地 影 模 
型 的 预报 结果 差异 逐渐 增 大 . 位 置 误差 在 一 天 内 可 
以 积累 到 几 个 毫米 , 一 些 误差 较 大 的 轨道 可 以 积累 
到 1 cm 左右 . 我 们 注意 到 在 4 条 轨道 中 , LEO3 的 位 
置 偏 差 显 著 地 大 于 其 他 3 条 轨道 , 这 是 因为 LEO3 轨 
道 进出 地 影 次 数 较 多 , 所 以 地 影 模 型 对 它 的 影响 更 


显著 . 


40 50 60 70 


Time/h 


图 11 不 同 地 影 模型 对 4 条 轨道 预报 精度 的 影响 


Fig.11 The influence of Earth shadow model on the accuracy of the four orbits 


3.5 ”潮汐 形变 报 动 、 后 牛顿 效应 和 大 行星 引力 摄 
动 对 轨道 预报 精度 的 影响 

潮汐 形变 摄 动 、 后 牛顿 效应 和 大 行星 引力 摄 

动量 级 相对 于 前 面 的 摄 动 都 比较 小 在 潮汐 形变 

摄 动 中 , 考虑 太阳 和 月 球 使 得 地 球 发 生 的 固体 潮汐 

S 变 ， 且 考虑 到 三 阶 近似 .相对 于 日 月 引力 , 太 
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阳 系 内 的 其 他 行星 的 引力 摄 动 较 小 , 另外 , 当 行 星 
与 地 球 的 相对 位 置 变 化 时 , 其 摄 动 效果 也 会 发 生变 
化 , 因此 , 各 个 行星 的 引力 摄 动 的 量 级 随时 间 会 发 
生变 化 . 一 般 情况 下 , 比较 显著 的 是 木星 的 引力 , 除 
此 之 外 , 水星、 金星、 火星、 土星 的 引力 摄 动 比 木 
星 的 引力 摄 动 小 一 到 两 个 量 级 , 天 王 星 和 海王 星 的 
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引力 摄 动 比 木星 小 3 个 量 级 或 更 多 . 这 里 研究 大 行 
星 引 力 摄 动 对 轨道 预报 的 影响 时 , 我 们 对 比 的 是 动 
力学 模型 是 否 包括 木星 、 水 星 、 人 金星 、 火 星 、 土 
星 的 引力 摄 动 的 预报 结果 . 关于 后 牛顿 效应 , 我 们 
考虑 的 是 地 球 引 力 场 的 广义 相对 论 效 应 修正 . 由 于 
广义 相对 论 下 的 动力 学 与 牛顿 体系 下 的 动力 学 差 
别 较 小 , 所 以 考虑 后 牛顿 效应 相当 于 再 增加 一 个 摄 
动 源 , 我 们 这 里 只 考虑 参数 化 的 后 牛顿 效应 的 一 阶 
项 HU. 我 们 计算 的 轨道 包含 2.2 节 所 述 的 各 个 类 型 ， 
总 数 2000 条 . 

图 12 展 示 了 潮汐 形变 摄 动 、 后 牛顿 效应 和 大 
行星 引力 摄 动 对 轨道 预报 精度 的 影响 , 图 中 每 个 点 
代表 一 个 空间 目标 , 方形 散 点 和 圆 形 散 点 分 别 代 表 
近 圆 轨道 和 大 偏心 率 轨道 , 纵 坐 标 表示 一 天 内 的 位 
置 偏差 . 作为 由 地 球 产 生 的 摄 动 , 潮汐 力作 用 和 后 
牛顿 效应 的 摄 动 作用 都 是 随 轨道 的 增加 而 减 小 的 ， 
而 作为 来 源 于 地 外 的 摄 动 , 大 行星 引力 摄 动 效果 随 
轨道 的 升 高 而 增 大 . 潮汐 形变 摄 动 在 一 天 内 可 以 使 
低 轨 空 间 碎 片 积累 百 米 以 内 的 位 置 误 差 , 而 对 于 地 
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球 同步 轨道 , 位 置 误差 在 1 m 以 下 . 所 以 对 于 地 球 同 
步 轨道 的 米 级 精度 和 中 低 轨 的 百 米 级 精度 , 潮汐 形 
变 摄 动 是 不 可 或 缺 的 . 后 牛顿 效应 积累 的 误差 比 测 
汐 形 变 摄 动 小 一 到 两 个 量 级 左右 , 对 于 中 低 轨 空间 
碎片 , 一 天 会 积累 10 m 以 下 的 位 置 误差 , 对 于 地 球 
同步 轨道 , 误差 则 减 小 到 0.1 m 左 右 . 因此 对 于 中 低 
轨 空 间 碎片 轨道 预报 的 米 级 精度 以 及 地 球 同 步 轨 
道 的 分 米 级 精度 , 它 是 需要 考虑 的 . 大 行星 引力 摄 
动 在 一 天 内 对 空间 碎片 轨道 预报 的 影响 在 0.1 m 以 
内 , 对 于 低 轨 空间 碎片 甚至 在 1 mm 以 内 , 它 的 摄 动 
效果 量 级 较 小 , 只 有 在 精度 要 求 十 分 高 的 场景 中 才 
需要 考虑 . 值得 注意 的 是 , 无 论 对 于 哪 一 种 摄 动 , 中 
轨 大 偏心 率 轨道 的 空间 碎片 的 误差 都 较 大 而 且 在 
图 中 分 布 较为 分 散 , 这 是 因为 当空 间 碎片 的 偏心 率 
较 大 时 , 它 所 处 的 轨道 高 度 分 布 也 会 比较 广泛 . 以 
闪电 轨道 为 例 , 它 的 偏心 率 达 到 0.74,， 这 使 得 它 的 
远地点 比 地 球 同 步 轨道 更 高 , 而 它 的 近地点 又 处 于 
低地 球 轨道 . 
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辐 12 ”潮汐 形 变 ( 左 上 )、 后 牛顿 效应 (右上 )、 大 行星 引力 (左下 ) 对 轨道 预报 精度 的 影响 


Fig.12 Effects of tidal deformation (upper left), post-Newtonian effect (upper right), and gravitational pull of large planets 


(lower left) on orbit prediction accuracy 
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动力 学 模型 选择 的 建议 
动力 学 模型 的 选择 是 在 预报 精度 与 预报 速度 
之 间 的 取舍 . 详尽 、 复 条 的 动力 学 模型 可 以 使 得 轨 
道 预报 的 精度 十 分 高 , 但 是 当面 临 大 批量 的 空间 碎 
片 时 , 预报 速度 是 同样 重要 的 需求 . 从 图 3 与 图 4 中 
我 们 已 经 看 到 , 同样 精度 需求 下 , 不 同 轨道 高 度 的 
空间 碎片 适用 的 引力 场 模型 的 阶 次 有 着 十 倍 甚至 
更 高 的 差别 , 而 阶 次 差别 在 十 倍 以 上 意味 着 在 引力 
场 的 计算 中 的 计算 量 差别 在 百倍 以 上 . 同样 , 日 月 
历 表 的 选择 、 潮 汐 形 变 的 考虑 与 否 等 都 对 轨道 预 
报 的 计算 量 有 着 重要 影响 

在 本 节 中 , 我们 会 考虑 各 种 轨道 预报 场景 ( 主 
要 是 不 同 的 轨道 类 型 和 不 同 的 精度 需求 ), 给 出 最 
优化 的 动力 学 模型 建议 .建议 的 动力 学 模型 可 以 
满足 精度 需求 , 同时 , 它 是 最 简化 的 、 计 算 量 需 求 
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已 起 取 间 


最 小 的 动力 学 模型 . 如 图 13 所 示 是 动力 学 模型 建 
议 图 , 建议 图 共 分 成 6 个 部 分 , 分 别 表 示 6 个 不 同 的 
精度 需求 , 从 最 高 的 1 cm 级 的 精度 到 最 低 的 1 km 级 
的 精度 . 纵 坐 标 为 不 同 的 轨道 类 型 , 除了 在 2.2 节 中 
已 经 说 明 的 将 空间 碎片 的 轨道 类 型 大 致 分 为 4 个 类 
型 外 , 对 于 每 一 个 类 型 ， 又 通过 高 度 的 差别 分 出 了 
更 精细 的 类 别 . 其 中 LEO 按 照 高 度 分 出 了 4 个 层级 ， 
MEO 按 照 高 度 分 出 了 3 个 层级 , 这 样 做 的 原因 是 即 
使 同样 归 为 LEO 的 轨道 , 使 用 的 动力 学 模型 依然 有 
很 大 差别 , 例如 在 引力 场 的 精细 程度 上 ， 几 乎 高 度 
每 提高 100 km 左右 ， 阶 次 就 可 以 有 一 定 的 下 降 ( 详 
见 图 3 与 图 4). 横 坐 标 列举 了 本 文中 研究 的 6 种 摄 动 
力 , 它们 分 别 是 地 球 的 非 球 形 引 力 位 、 大 气 阻 力 、 
日 月 第 三 体 引力 摄 动 、 光 压 摄 动 、 潮 汐 形 变 摄 动 、 
后 牛顿 效应 和 木星 等 天 体 的 第 三 体 引 力 摄 动 . 图 中 
的 引力 场 用 数字 表示 , 例如 数字 10 即 表示 引力 场 采 
910 x 10 的 阶 次 , Y (绿色 背景 ) 表 示 考 虑 此 项 摄 动 ， 
而 N (红色 背景 ) 表 示 不 需要 考虑 此 项 摄 动 . 另外 , 在 
行星 历 表 ( 日 月 第 三 体 引力 摄 动 ) 一 列 中 Ana (绿色 
背景 ) 表 示 使 用 简化 的 分 析 历 表 ， 分 析 历 表 的 内 容 


~ 


月 所 
已 在 前 文中 说 明 . Num (红色 背景 ) 表 示 需 要 使 用 数 
值 历 表 , 对 JPL 历 表 来 说 ， 只 要 使 用 其 不 太 老 旧 的 
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版 本 即 可 (比如 DE4** 系 列 都 没有 问题 ). 动力 学 模 
型 建议 图 中 并 没有 涉及 地 影 模 型 ， 因为 从 3.4 节 来 
看 , 我 们 设 定 的 所 有 场景 中 , 柱 形 地 影 近似 造成 的 
误差 都 是 可 以 忽略 的 . 另外 , 由 于 地 影 模 型 的 差异 
在 步 长 为 1 s 左 右 或 更 短 的 轨道 预报 中 才 会 有 所 体 
现 , 而 空间 碎片 的 轨道 预报 通常 不 会 采取 如 此 小 的 
步 长 , 所 以 柱 形 地 影 带 来 的 误差 不 会 对 空间 碎片 的 
轨道 预报 造成 影响 . 

举例 来 说 ， 若 有 一 处 于 900 km 高 度 的 低 轨 空 
间 碎 片 ， 我 们 需要 对 它 进 行 一 天 内 位 置 误差 小 于 
1 km 的 轨道 预报 , 那么 根据 图 13, 推荐 的 动力 学 模 
型 应 当 是 考虑 3 种 摄 动因 素 , 分 别 是 15x15 阶 次 的 非 
球形 引力 位 、 大 气 阻 力 、 分 析 历 表 下 的 日 月 第 三 体 
引力 . 若 有 一 处 于 地 球 同步 轨道 的 空间 碎片 , 我 们 
需要 对 它 进行 一 天 内 位 置 误差 小 于 10 m 的 轨道 预 
报 , 则 推荐 的 动力 学 模型 是 考虑 3 种 摄 动因 素 , 分 别 
是 5x5 阶 次 的 非 球形 引力 位 、 数 值 历 表 下 的 日 月 第 
三 体 引 力 和 光 压 . 

在 动力 学 模型 建议 图 中 , 轨道 类 型 的 划分 依据 
九 道 高 度 与 轨道 偏心 率 , 而 不 是 卫星 的 功能 
源 , 因为 轨道 半 长 径 和 轨道 偏心 率 是 决定 适用 
动力 学 模型 的 主要 因素 . 对 于 某 一 类 型 的 轨道 ， 
我 们 可 以 通过 它们 的 轨道 高 度 与 偏心 率 判 断 它 们 
处 在 动力 学 模型 建议 图 的 哪 一 部 分 . 例如 太阳 同步 
轨道 是 一 种 主要 处 在 低 轨 的 小 偏心 率 轨道 , 地 球 扁 
率 摄 动 使 得 其 升 交 点 经 度 的 进 动 速度 与 地 球 的 公 
转速 度 相 同 ,在 空间 上 表现 为 其 轨道 高 度 与 轨道 
倾角 满足 一 定 关系 . 不 同 的 太阳 同步 轨道 所 处 的 
轨道 高 度 不 同 , 那么 建议 的 动力 学 模型 也 不 同 , 比 
如 中 国 第 1 颗 太 阳 探 测 科 学 技术 试验 卫星 “ 尺 和 号 ” 
处 在 一 条 轨道 高 度 约 为 520 km 的 太阳 同步 轨道 上 
如 果 我 们 希望 对 它 进 行 轨道 预报 的 精度 控制 在 位 
置 误 差 每 天 不 超过 1 km, 那么 本 文 推荐 的 动力 学 模 
型 如 图 13 第 1 行 中 间 所 示 . 再 例如 , 截至 2022 年 5 月 ， 
星 链 计划 的 所 有 在 轨 卫 星 均 处 在 高 度 500 km 到 
600 km 之 间 ， 那 么 它们 都 可 以 被 归 类 到 图 13 的 
第 1 行 中 , 即 轨 道 高 度 在 300 km 到 600 km 的 LEO 轨 
道 , 遵循 相应 的 动力 学 模型 建议 
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13 动力 学 模型 建议 图 . 图 中 6 个 子 图 分 别 表示 6 种 精度 需求 , 每 一 行 对 应 一 个 轨道 高 度 区 间 , 高 度 区 间 在 左 侧 的 表格 中 说 明 . 图 中 的 Y 表 示 动 力学 
模型 应 考虑 该 摄 动 , N 表 示 动 力学 模型 不 需 考虑 该 摄 动 . 引力 场 (Gravitational field) 这 列 的 数字 表示 应 采用 的 引力 场 阶 次 . 


Fig.13 Suggestions on the reference dynamic models. The six subplots in the figure represent six different accuracy 
requirements. Each row corresponds to an orbital altitude range, which is explained in the table on the left. In the figure, “Y” 
indicates that the perturbation should be considered in the dynamic model, while “N’” indicates that it does not need to be 


considered. The numbers in the “Gravitational field” column indicate the order of the gravitational field that should be used. 


5 总结 和 展望 一 个 参考 . 轨道 预报 的 结果 表明 , 不 同 摄 动力 的 量 
在 本 文中 , 我 们 使 用 基于 RKF7 (8) 的 数值 法 轨 ”级 差距 民 大 , 同时 , 同一 担 动 力 在 不 同 的 预报 场景 


中 


道 预 报 计算 了 大 批量 的 空间 碎片 轨道 , 并 且 通 过 使 《中 发 挥 的 作用 也 有 巨大 差别 . 挑选 正确 的 动力 学 模 


;办 


] 不 同 的 动力 学 模型 来 评估 各 个 摄 动力 对 轨道 预 ”型 将 不 仅 会 使 轨道 预报 速度 有 质 的 飞跃 , 而 且 可 以 
报 精度 的 影响 . 在 以 往 的 实践 中 , 我 们 仅仅 知道 某 ”减少 轨道 预报 过 程 对 外 部 的 依赖 

项 摄 动力 的 量 级 大 概 在 什么 范围 , 而 并 没有 普 适 空间 碎片 的 数目 巨大 、 轨 道 类 型 各 异 , 预报 不 
和 具体 的 数据 , 本 文中 的 计算 结果 可 以 为 大 家 提供 同 的 空间 碎片 需要 使 用 不 同 的 策略 , 而 传统 的 轨道 
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预报 是 需要 人 为 选择 动力 学 模型 的 . 如 果 对 于 每 一 
个 或 每 一 种 空间 碎片 人 工 选择 动力 学 模型 , 这 将 是 
十 分 不 方便 的 , 而 如 果 粗 糙 地 选择 同一 种 动力 学 模 
型 ,精度 和 预报 效率 又 难以 得 到 保障 . 我 们 的 结果 
也 许 能 为 这 个 问题 提供 一 个 可 行 的 解决 方案 , 基于 
本 文 的 结果 , 可 以 为 不 同 的 轨道 类 型 、 不 同 精度 要 
求 的 空间 碎片 自动 分 配 最 优 的 轨道 预报 方案 , 可 以 
极 大 地 节省 人 力 和 算 力 的 成 本 . 

动力 学 模型 建议 图 既 可 以 作为 轨道 预报 时 动 
力学 模型 选择 的 参考 , 也 可 以 用 于 特定 动力 学 模型 
下 的 预报 误差 的 估计 . 我 们 尽 可 能 地 在 有 限 的 图 表 
中 列 出 尽 可 能 多 的 预报 场景 , 但 是 真实 的 预报 场景 
会 更 加 具体 和 复杂 , 我 们 会 在 后 续 使 预报 场景 更 加 
细 化 , 并 且 提 供 更 加 方便 的 使 用 方法 , 例如 选择 动 
力学 模型 的 源 代码 , 自动 选择 动力 学 模型 的 轨道 预 
报 程序 等 . 
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Reference Dynamic Models for Space Debris Orbit Prediction 
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ABsTRAcT In different orbit prediction scenarios, the dynamic models used are also different. For exam- 
ple, the atmospheric drag is a very important perturbation force in the prediction of space debris in low 
orbits, but at mid-to-high orbits, the atmospheric drag can be ignored. In order to select the most suitable 
dynamic models for orbit prediction in different scenarios, this article provides a large number of orbits 
for different accuracy requirements and different orbit types. By comparing the prediction results with 
different dynamic models, the optimal dynamic models for various prediction scenarios are recommended. 
This work can serve as references and standards for the selection of reference dynamics models for different 


orbital types of space debris in orbit prediction. 
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